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Abstract 



We search for rational, four-dimensional maps of standard type (x n+ i — 2x n + 
x n -x — ef(x, e)) possessing one or two polynomial integrals. There are no non- 
trivial maps corresponding to cubic oscillators, but we find a four-parameter family 
of such maps corresponding to quartic oscillators. This seems to be the only such 
example. 



Suris [1] has found all one-degree-of-frccdom Hamiltonian systems that possess inte- 
grable discretizations of standard type. If f(x,0) is rational, it must be a polynomial of 
degree < 3 and f(x,e) is rational. Naive counting suggests that if some 2-d.o.f. continu- 
ous systems possess a second integral, then some 2-d.o.f. symplcctic maps might possess 
one integral. Here we search for such maps, which turn out to be rare: we find one four- 
parameter family of maps with one integral, a three-parameter subset of which is rotationally 
invariant and hence has a second integral; this subset corresponds to discretizations of the 
two-parameter family of Hamiltonians \{p\ +P2) + A(x\ + x\) + B{x\ + x 2 ) 2 . The rarity is 
perhaps because we want a whole family of integrable maps (depending on e), not just an 
isolated map. 

A symplectic map of standard type is written 

x„+i - 2x„ +x„_i = ef(x„,e), f (x, e) = W(x, e) (1) 

and the small parameter e may be thought of as the square of the time-step in a discretization 
of x = f (x, 0). Assume that f and V have expansions 

00 00 

f (x, e) = j*i W, V(x, e) = E ejV o W 
j=o j=o 

We want an integral <I>(x, y, e) such that 

$(x n+ i,x n ,£) = $(x n ,x n _i,£) for all n (2) 

and assume that 

$(y,x,e) = $(x,y,e) = $o(x,y) + e$i(x,y), 

the first restriction following from the reversibility (under n ^ —n) of (1). 

Our construction follows that of Suris. Define u„ = x„ — x„_i +ef (x„, e)/2; then from 
(1) and (2) the condition for $ to be an integral becomes 

$(x,x + ef(x,e)/2 - u,e) = $(x, x + ef (x, e)/2 + u, e) (3) 

which may be expanded as a Taylor series in e. At order e°, we get $o( x :y) = ¥'( x — y) 
(where ip is necessarily even), and at order e, 

(V^(u)) • f (x) = *i(x,x-u)-$i(x,x + u). (4) 

In the one-degree-of-freedom case, Suris's method for solving the functional equation (3) is 
to: 

(i) Differentiate the order e term twice, the order e 2 term once, and combine with the order 
e 3 term; this gives six equations in the six unknowns 3>^(x, x ± u) (i = 1, 2, 3) which 
turn out to have rank 5. The consistency condition for these equations is separable 
in x and u, and he obtains <p"' (u) / 'tp' (u) = c, giving three candidates for ip: u 2 /2, 
(1 — cos oju) I 'lo 2 , and (cosham — 1)/oj 2 - 

(ii) For each candidate, deduce the functional form of $1 from (4). 

(iii) Take the most general function of this form, substitute into (3), and solve for f(x,e). 
In each case, the candidate for <p did in fact give a family of solutions for /. With 

two degrees of freedom, the same steps could be followed: for (i), to get an overdetermined 
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system it is necessary to differentiate the order e term three times, etc., proceeding to the 
order e 4 term; this gives a rank- 18 set of twenty equations in twenty unknowns and two 
consistency conditions relating y)W(u) and f^ J '(x); these equations are complicated and 
probably not separable. However, the severest restriction is the new one that f(x, e) must 
be a gradient: 

(iv) Given $1, require |^ = 

This must be checked for any proposed (p. We have taken ip(u) = {u\ + u 2 )/2 and 
searched for solutions of (3) satisfying (iv). From (4), <I>i(x, y) must be a polynomial, and 
at most quadratic in each variable, corresponding to a quartic potential Vo- Hence we take 

2 

$i(x,y)= Pvkix\x{yly l 2 {Pijki=Pkiij Vi,j,k,l). 

i,j,kj=0 

The consistency conditions on f when solved directly from (3) are very complicated, so we 
determine them term-by-term. Once a solution at the first few orders is obtained we go 
back and check (3) directly. First determine Vo from (4): 

Vo = -(2x 2 poioo + x\ (poioi + 2p 20o) + 2x 2 1 po20i + x 2 p 202+ 

xi (2piooo + x 2 (2piooi + 2pnoo) + x\ (2pi 00 2 + 2pn i) + 2x3^1102) + 
x\ (pioio + 2 P2000 + x 2 (2pmo + 2p 2 ooi) + x\ (pun + 2^2002)) + 

x\ (2?52010 + 2x 2 P201l) + ^P2020) 

and then at order e n (n = 2, 3, 4) we have 

u.(f„_i + + h n (x))=0 (5) 

where g„(x) is a polynomial and g 2 = 0. Solving (5) for f„_i gives a consistency condition 
at that order and a value which can be substituted into the next order. The consistency 
conditions take the form of polynomials in x\, x 2 (whose coefficients are functions of the 
Pijki) which must be zero; this gives polynomial equations, of degree n, in the p's. To simplify 
the equations we took various choices for Vo and attempted to solve these consistency 
equations. (The actual equations and solutions are derived using a symbolic manipulator.) 
Without any other restrictions there are 14 equations at order 2, 27 at order 3, and 38 at 
order 4. 

If the map is to be completely intcgrable, then the vector field it approximates as 
e — > must be also; hence we first tried Vb's corresponding to known integrable oscillators 
[2,3]: Ax\ + Bx 2 - x\x 2 - 2x 2 ; A(x\ + 16x 2 ) - x\x 2 — I6X2/3 (two cases of the Henon-Heiles 
Hamiltonian); x\x 2 + 2x\\ x\ J r l2x\x\- s r l&x 2 ; x\ + 6x\x% + 8x 2 ; and A{x\ + x 2 ) + (x\ + x 2 ) 2 . 
Only in the last case is there a consistent solution for the p's (given below), although 
sometimes it is necessary to continue to order e 4 to reach the inconsistency. 

Now consider whether there are any maps with only one integral. If so, then (4) 
shows (take d/du and set u = 0) that the integral must correspond to the Hamiltonian as 
e — > 0, i.e. $i(x,x) = — Vo(x). We are not interested in the trivial solutions — those which 
give linear maps, those which give uncoupled maps, or those which uncouple after a linear 
symplectic change of variables. For the form 

Vo = Ax\ + Bx\ + ax\ + (5x\x 2 + ^x x x\ + 8x1 
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the only such solutions found had A = B and 3(aj + pS) = [3 2 + j 2 ; a long calculation shows 
that these maps all uncouple after a linear symplectic change of variables. Thus there do 
not seem to be any nontrivial maps of this form with one integral corresponding to cubic 
oscillators. 

For the form 

V = Ax\ + Bx\ + ax\x 2 + (ix\x\ 

we found no nontrivial solutions. 
For the form 

Vb = Ax\ + Bx\ + ax\ + (lx\x\ + -fx 2 

we found a four-parameter family of nontrivial solutions with (i = 2a = 2-f and arbitrary 
A, B. This seems to be the only nontrivial solution: Let 

Q = b ((1 - ce)x\ + (1 - de)x 2 2 ) . 

Then 



V = ln(2(l - ce)(l - de) -sQ)- ^- 



be 2 v v ' y 1 ^' e 
(a + 2d)(l - ce) + Q , (a + 2c) (1 - de) + Q 

Jl=Xl— 3—7- r — — , J2=X 2 - 



(6) 



{l-ds){l-ce)-eQ/2' J (1 - <fe)(l - ce) - sQ/2 

with integral 

-2$j = a(x iyi + x 2 y 2 ) + b(x lVl + x 2 y 2 ) 2 /2 + c(x\ + y%) + d(x\ + y\) 

When c = d, the potential V is rotationally invariant and we expect a second integral; 
the method determines this automatically, for in this case one finds a term e{x\y 2 — x 2 y{) 2 
in $1 with V independent of e. The obvious extension of this map to n degrees of freedom 
is also integrable: 

_ a + 2c + b|x| 2 
_ l-ce-efe|x| 2 /2 X 

y = -^ln(2-fe|x| 2 )-^ ( 7 ) 
be z e 

, r a + 2c 2 b 4 

y = ^|x| +-|x| 

Grammaticos et al. [4] have proposed an integrability test for maps: movable singu- 
larities should not propagate in time, and memory of the initial condition should survive 
the singularity. To apply the test (for simplicity, for 2 degrees of freedom), write the map 
in the form x„ + i = — x n -i + ct/(l — /3 2 |x„| 2 )x„ and take initial conditions Xo arbitrary, 
xi = (a;ii,0) (this suffices because we can rotate coordinates). Compute the next three 
iterates and then let in — > this gives the iterates 

(x 00 ,x 01 ), (1//3.0), (00, -xoi), {-1/(3,0), (-zoo, zoi) 

showing that the proposed integrability test is satisfied here. 

We can now use the new "angular momentum" integrals to reduce the map to one 
degree of freedom. This illustrates that the two processes, (i) forming an integrable map 
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approximating a continuous system and (ii) reducing by the rotational symmetry, do not 
commute. Rotate coordinates so that Xi = yt = for i > 2. Write the map as 

x' = x + p 
p'=p + x' ff (|x'|) 

and define new variables (chosen to correspond to the continuous case) 

L =xip 2 - x 2 pi 
(r,9) : polar coordinates for (xi,x 2 ) 
p r = (X!P! + x 2 p 2 )/r 
Pe=L 

in which the map can be written (after much algebra) 

(r 2 )' = ^ + (p r +r) 2 
L 2 



( r Pr)' = ^I+ r Pr+Pl + (r') 2 9(r') 
= (r') 2 (l + g(r / ))-r(r+p r ) 
L' = L 

6' = (9 + sin- 1 A 



(8) 



with integral 



L 2 ( ( L 2 

r(r+p r )(a + br(r+p r )/2)+c[^- + (r + p r ) 2 +r : 



r-2 



However, starting with the continuous system corresponding to (7) (x = VVo(l x D) an d 
performing the corresponding reduction gives 

which does not have an integral discretization of standard type. So we have extended Suris's 
list of systems possessing integrable discretizations by one; the catch is the our reduced 
system (8) is not of standard type. This suggests that there may be rational integrable 
discretizations of any rational Hamiltonian if one enlarges the allowed class of maps. 
My thanks to James Meiss for pointing out [1] and for many helpful discussions. 
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